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MULTIPATH  VLF  PROPAGATION  EFFECTS  ON 
CORRELATION  RECEIVERS 


INTRODUCTION 

The  propagation  of  very-low  frequency  (VLF)  waves  over  terrestrial  paths  is  often  treated  by  con¬ 
sidering  the  space  between  the  earth  and  the  ionosphere  to  be  the  interior  of  a  waveguide,  with  the 
earth  as  a  lower  boundary  and  the  D  and  E  layers  of  the  ionosphere  as  an  upper  boundary.  Ordinarily 
several  waveguide  modes  can  propagate  with  low  attenuation  in  the  guide.  Where  two  modes  have 
approximately  equal  strength  and  arrive  with  approximately  180s  phase  shift  between  their  wave  com¬ 
ponents,  a  modal  interference  null  is  said  to  occur.  Typical  examples  of  such  modal  interference  nulls 
are  shown  in  Fig.  1,  taken  from  Rhoads  and  Gamer  [1],  The  predicted  curves  for  Fig.  1  were  those 
generated  using  the  theoretical  waveguide  formulation  of  Wait  and  Spies  [2], 


Fig.  I  —  Comparison  of  theoretical  results  from  Wait  and  Spies 
(21  with  observed  records  in  flight  from  San  Francisco  to  Hono¬ 
lulu,  May  18-19,  1965,  for  radiated  power  P,  -  1  kW 


When  an  airborne  facility  transmits  VLF  waves  using  a  trailing-wire  antenna,  a  complicated  time- 
varying  field  is  produced  at  a  receiver  [3,4].  The  inclined  antenna  can  be  considered  a  superposition  of 
an  elevated  vertical  and  horizontal  antenna.  Both  quasi-transverse  magnetic  and  electric  modes  are  gen¬ 
erated  by  the  inclined  antenna.  The  amplitudes  of  theae  modes  vary  with  the  angle  between  the  hor¬ 
izontal  component  of  the  antenna  and  the  direction  to  the  receiver.  If  the  horizontally  and  vertically 
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generated  resultant  fields  are  nearly  equal  in  magnitude  at  a  receiver  location,  rather  deep  oscillations  in 
the  received  field  strength  may  occur  during  the  period  of  an  orbit.  These  oscillations  often  occur  dur¬ 
ing  nighttime  propagation  conditions  but  can  occur  during  daytime  also.  Figure  2  shows  the  signal 
reception  pattern  versus  distance  and  time  for  an  airborne  transmitting  facility  flying  in  good  orbit  (ver- 
ticality  =70%)  during  daytime.  Figure  3  shows  the  same  pattern  for  a  poor  orbit  (verticality  =10%) 
during  nighttime.  The  observed  fluctuations  range  from  mild  to  extreme. 


Fig.  2  -  Field  strength  versus  distance  measured  during  day-  Fig  3  -  Field  strength  versus  distance  measured  during  nigh- 

light  hours  on  an  aircraft  flying  radially  toward  an  airborne  time  hours  on  aircraft  flying  radially  toward  airborne  broad¬ 
broadcasting  facility  in  good  orbit  near  Jacksonville,  Florida  casting  facility  in  poor  orbit  near  Jacksonville,  Florida  on 

on  March  27,  1971,  with  the  frequency  being  19.4  kHz  April  5,  1971,  with  the  frequency  being  25.6  kHz 

An  important  task  of  airborne  VLF  facilities  is  to  broadcast  wideband  spread-sprectrum  VLF  com¬ 
munication  signals.  These  signals  may  be  distorted  by  the  propagation  through  the  waveguide  and  by 
the  variation  in  aircraft  direction  in  orbit.  To  determine  the  deleterious  results  caused  by  these  cir¬ 
cumstances,  it  is  necessary  to  examine  the  correlation  of  the  received  signal  with  an  exact  replica  of  the 
transmitted  waveform  to  estimate  the  effective  signal  reduction. 

A  study  of  correlation  receiver  distortion  due  to  propagation  effects  was  performed  by  Rothmuller 
[SI  for  the  case  of  broadcast  from  a  fixed  transmitter  to  a  fixed  receiver.  Rothmuller  showed  that  when 
two  multipath  signals  interfere  at  a  point,  two  correlation  peaks  occur  separated  in  receiver  synchroniza¬ 
tion  time  r*  by  approximately  the  duration  of  a  spread-spectrum  signaling  element  (chip).  The  max¬ 
imum  value  of  the  correlation  could  be  on  the  earlier  or  later  correlation  peak.  In  the  present  report  we 
go  a  step  further  by  including  the  effects  of  aircraft  orbit  on  correlation  receiver  signal  output.  We 
study  the  condition  of  a  time-varying  wideband  multipath  situation  common  in  airborne  facility  usage 
by  considering  the  case  in  which  the  signaling  elements  (chips)  are  minimum-shift-keying  (MSK) 
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waveforms  and  the  communication  decision  is  based  on  the  coherent  integration  of  a  series  of  such 
waveforms. 

In  some  circumstances  the  performance  of  the  correlators  depends  greatly  on  which  instant  within 
a  chip  the  correlator  is  synchronized  to.  It  appears  probable  that  a  bit  decision  based  on  several  correla¬ 
tions  within  a  chip  period  would  be  superior  to  a  decision  based  on  a  single  correlation  within  a  chip. 

THEORY 

Elementary  Signal  Theory 

Figure  4  shows  schematically  the  outlines  of  the  communication  system  under  discussion.  Figure 
S  shows  a  conceptual  version  of  a  receiver  and  demodulation  portion  of  Fig.  4.  The  operation  of  the 
integrator  portion  of  the  equipment  in  Fig.  S  is  to  sample  the  received  waveforms  A  (t)  and  to  multiply 
them  by  the  desired  waveform  A( (i)  and  A2(/)  and  to  integrate  over  a  time  equal  to  a  bit-keying- 
element  period  TB,  centered  about  the  time  t,  the  midpoint  of  the  /  th  expected  bit  interval.  The 
transmitted  waveform  is  constructed  as  a  linear  superposition  of  A(0)  and  A2(r)  according  to 

h(t)  -  Xt  A | (/)  +  Y,  h2(t).  (1) 

The  outputs  of  the  integrators  will  yield  the  values  of  X,  and  Y,  which  contain  the  transmitted  informa¬ 
tion, 

x,  -  (Ta)-‘  S^rii  dt  (2) 

and 

((^Vj  h2U)hU)  dt,  (3) 

as  long  as  AiO)  and  h2U)  are  orthonormal  to  each  other: 

i  C  ^ g/2 

( Tb)~ '  yrTm  M»>  M'>  *  -  *</.  (4) 

Here  8y  is  the  Kronecker  delta  function,  equal  to  zero  for  i  ^  J  and  equal  to  1  for  /  -  j. 

Correlation  of  Orthogonal  Waveforms 

A  set  of  orthogonal  wideband  waveforms  A,(r)  and  h2(t)  can  be  generated  by  using  a  baseband 
waveform  hB(t)  to  modulate  both  a  sinusoidal  and  cosinusoidal  wave: 


A^r)  -  hBU)  cos  <o0t. 

(5) 

h2(t)  -  As(r)  sin  w0t. 

(6) 

hB(t)(Xi  cos  w0t  +  Y,  sin  & i0t). 

(7) 

NOISE 


VIE  PROPAGATION  CHANNEL 

Fif .  4  —  Model  of  the  communication  ayitem 
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Fig.  5  —  Model  of  the  receiver-demodulator 


The  values  of  X/  and  Yt  give  the  strength  of  the  cosinusoidal  and  sinusoidal  components.  A  correlation 
receiver  can  then  be  made  to  extract  the  value  of  X ,  and  Y,: 

Xi  "  2(ra)  J tr  t g/i  hUHXi  cos«o /  +  Y,  sin  w0t)  hB(t)  cos  w0t  dt  (8) 

and 

r't+V2 

“  2(7fl)  1  J,_T^2  M*)  (X,  cos  a »0f  +  Y,  sin  u0l)  hB(t)  sin  w0t  dt.  (9) 

A  propagation  delay  time  of  tp  modifies  the  correlator  outputs  to  give  new  ones  which  we  call  T,(t.) 
and  Y,(t „). 

r'i*  V2 

X,(tp)  -  2(Tb)~x  JlrTi/1  hB(t  +  rp)  (*,  cos  <o0(f  +  tp)  +  Y,  sin  «u0(f  +  t,)1 

x  hB(t)  cos  a i0t  dt  (10) 

and 

/.<,+  ig/i 

irTgii  +  cos  ot0u  +  rp)  +  Y,  sin  w0(t  +  t,)] 

x  hB(t)  sin  w0t  dt.  (U) 

We  find  that  the  vector  fl/(rp)  —  ( X,(tp ).  K,(r.))  is  approximately  the  phase-routed  version  of 
R,  -  ( X/.Y ,).  We  have 

Xib,)-  Ta  (X,  cos<u0t„  +  Y,  sincuoT,)  J  hB(t )  hB(t  +  rp)  dt  (12) 


/•(l+V2 

Y,hp)  -  Tg  l(~Xi  sin  o>0tp  +  Y,  cos  warp)  j,rTg/J  hB(t)hBU  +  rp)  dt.. 
Thus  we  can  write  this  equation  in  vector  form  if  Rt  —  (Xt,  Y,): 


R,(rp)  -  TB 


,[  (COS 

|XrV2  f>B(t)hg{t  +  rp)d^_  sil 


os  w0rp  sin  cuqt  p  _ 
sin  coqt p  cos  woTp  R'' 


The  magnitude  of  the  vector  Rt( rp)  —  [Xi(rp)2  +  K,(t#)21i/j  is  given  by 

l*,(r,)l  “  Tj'\K,\  hBU)hBU+rp)dt.  (15 

The  phase  of  the  vector  Ri(tp )  is  given  by 

$/(»>)  -  arg(Xiirp)  +  iY,(rp))  -  arg(X,  +  iYt)  +  war p  (16 

Thus,  the  magnitude  of  the  R/(rp)  vector  varies  slowly  with  rp,  and  the  phase  <h(rp)  changes  rapidly. 
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Correlation  of  Other  Orthogonal  Waveforms 

A  more  complicated  set  of  orthogonal  waveforms  can  be  generated  using  two  orthogonal  baseband 
waveforms  hBl(t)  and  hB2(t)  such  that  for  /,  j  -  1,2 

(Tb)-'  hmit)hBjit)  dt  -  8 0,  (17) 

where  htJ  is  the  Kronecker  delta  function.  The  transmitted  waveform  hit)  is  given  by 

hit)  -  X,hB\it)  coso»0r  +  Y,hB2it)  sin  (18) 

To  extract  the  value  of  X,  and  Y,  from  this  waveform,  we  need  a  correlator  receiver  to  do  the  following 
integrations: 

/•'/+ VJ 

T(-  2(ra)-'  j,  lXthBiit)  cos  w0l  +  Y,hB2it)  sin  <o0t][/ifl,(t)+/is2(t)I  coso>0f  dt  (19) 
and 

r  'i*  V2 

Y,~  2iTB)  1  J  jyj  [A'/AAi<r )  coscu0t  +  Y,hB2it)  sin o*0/l [A0i (r)H-A93(f )]  sin<u0r  dt.  (20) 
When  the  time  of  arrival  is  changed  to  /  —  t  +  rp,  we  have 

rii+TJ2 

X/ir,)  -  2(TB)  1  J,rTji  \XihB\it  +  rf)  coso»0(f  +  rp)  +  Y,hB2it  +  tp)  sin«0(»+Tp)] 

x  lhj|(/)  +  hB2lt)]  cos  tug t  dt  (21) 

and 

r'/+VJ 

”  2iTB)  1  JlrT^2  lXihBi(t  +  rp)  coso) 0(f  +  Tp)  +  Y,hB2it  +  t p)  sin  <o0(t  +  t,)] 


x  lhB,(t)  +  hB2(t)l  sin  o>0t  dt. 


This  gives 


X,(tp)  -  X,  cos  w 0tpTb  1  J  hB,(t  +  Tp)[hai(r)  +  hB2(t)]  dt 
+  Y,  sin  «oTp Tb  1  J" hB2(t  +  rp)\hBXit)  +  hB2(t) )  dt 


,  r'/+V2 

Yii* p)  "  ~ Xt  sin  ioqt pTB  i2  hB\(t  +  r p)lhB\(t)  +  hB2it)]  dt 
rl,*Tg/2 

+  Y,  cos  wqt p Tb  J,  jy 2  hB2it  +  Tp)  I hB](t)  +  hB2ft))  dt. 
If  the  two  baseband  modulation  functions  are  sufficiently  orthogonal  that 

/•  r.+ TJ2 

((_V2  hBi(t+Tp)hB2(l)*  dt^O 

and  sufficiently  similar  that 

Tbx  Xr7y2  +  Tp)*ai h)  *  “  n'  Sif-Tf/2  ^ ^  +  hB20)  dt. 

we  again  obtain 

,  r'/+V2 

Xi(tp)  -  (X,  cos  «<fTp  +  Yi  sin  m<jrp)TB  1  f  hB\it  +  rp)hBiit)  dt 
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and 

ft.+  TJl 

y/(v  “  (-  X,sinat0Tp  +  Y,cosu0T,)Ti'  Jt-Tg/i  hiit+r  p)  hBX(t)  dt.  (28) 

When  the  more  complicated  correlators  of  Eqs.  (19)  and  (20)  are  used,  and  the  orthogonal  waveforms 
satisfy  the  conditions  of  Eqs.  (25)  and  (26),  the  type  of  waveform  of  Eq.  (18)  gives  a  behavior  similar 
to  that  of  Eq.  (7). 

Modulation  Methods 

The  values  of  Rt  can  be  modified  at  thejransmitter  to  initiate  the  communication  of  information. 
In  amplitude  modulation  the  magnitude  of  R,  is  varied.  In  phase-shift  keying  the  phase  angle  <£,  is 
modulated.  Thus,  in  general,  information  may  be  transmitted  by  changing  the  transmitted  vector  R,  in 
some  prescribed  way.  The  stability  of  the  received  vector  S,irp)  in  amplitude  and  phase  is  important  in 
the  performance  of  communication  systems.  In  the  realistic  communication  channel  the  values  of 
R/(rp)  will  also  depend  on  the  noise  present.  In  this  study  we  neglect  noise  and  simply  study  the 
effects  of  propagation  on  Riirp). 

MSK  Signals 

In  the  field  of  VLF  communications,  minimum  shift  keying  (MSK)  is  a  commonly  used  type  of 
spread-spectrum  modulation.  MSK  is  mathematically  described  in  terms  of  a  transmitted  waveform 
sit),  a  baseband  modulation  uit),  and  a  chip  waveform  pit),  as  follows: 


s(/)-Re[u(f)e/“0'], 

(29) 

N  . 

r  1 

«(/)-  £ 
*-0 

[*//„  Pit-  nT)  -  i  Y,  Q„  p  it-T/2  -  «7’)J. 

(30) 

pit)  -  V?[cos(t 

-  0.  1 1 1  >  T/2. 

This  is  a  specific  version  of  the  signal  described  in  general  by  Eq.  (18),  in  which  the  orthogonal 
baseband  waveforms  hBi it)  and  hB1it)  are  given  by 

T-1 

hBiit)  -  VT  £  lnpit  -  nT)  (32) 

»  “0 

and 

?" 

hBJit)  -VT  £  Q'PU-  T/2  -  nT).  (33) 

n-0 


-  T/2  <  /  ^  T/2. 


(31) 


In  these  above  equations  /„  and  Q„  are  pseudorandomly  chosen  integers  of  value  +1  or  -1.  It  is 
assumed  that  both  the  transmitter  and  receiver  possess  synchronized  pseudorandom-number-sequence 
generators  that  provide  the  I„  and  Q„  sequences.  The  frequency  tag  is  the  carrier  frequency  of  the 
broadcast.  The  time  interval  T  is  the  duration  of  a  chip.  Since  two  chips  are  transmitted  during  each 
chip  duration,  the  chip  length  re  is  T/2.  The  MSK  modulation  method  gives  rise  to  a  constant- 
amplitude  waveform  with  continuous  phase  as  long  as  X,  and  Y,  are  equal  to  ±  1 .  The  bandwidth  of 
the  transmitted  signal  is  inversely  related  to  the  chip  period.  The  value  of  the  waveform  sit)  depends 
on  one  value  of  /„  and  one  value  of  Q„  simultaneously  except  at  the  instant  when  a  transition  is  being 
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made  from  one  value  of  /,  or  Q„  to  the  next.  At  the  transition  time  the  amplitude  of  the  waveform 
undergoing  transition  is  zero,  because  the  pit)  function  goes  smoothly  to  zero  at  the  beginning  and  end 
of  its  range.  At  that  instant  the  value  of  sit)  depends  only  on  the  other  integer  Q„  or  /„.  Let  us  con¬ 
sider  the  autocorrelation  function  B(j )  for  the  MSK  waveform  s(r): 

Um  fj  sit) sit  +  r)  dt.  (34) 

Use  of  sit)  from  Eq.  (29)  will  make  B(r)  nonzero  as  long  as  |r|  <  T.  When  It  I  >  T,  the  contribu¬ 
tions  to  B(r)  will  all  be  from  terms  having  the  products  Q„  Im,  Qn  Qm,  and  /„  Im  with  n  m.  These 
terms  are  as  often  +1  as  they  are  -1  for  an  orthogonal  or  uncorrelated  set  of  {Q„  }  and  (/,)  ;  so  a 
summation  based  on  these  values  should  remain  small.  Thus  a  large  value  for  the  autocorrelation  func¬ 
tion  occurs  only  for  |r|  <  T,  when  Z,2  and  Q}  terms  contribute.  Since  1}  -  0„2  -  1 ,  these  terms  are 
easy  to  evaluate.  The  correlation  vector  is  closely  related  to  the  autocorrelation  function  B(t). 

At  very  low  frequencies  the  electromagnetic  waves  may  be  considered  to  propagate  by  distinct 
waveguide  modes  or  hops.  The  time  delays  between  dominant  modes  and  hops  are  small— less  than  a 
millisecond  [6]— producing  the  coherent  overlap  of  one  chip  with  its  time-delayed  image.  At  some 
receiver  locations  the  interference  of  the  two  waves  with  each  other  results  in  "selective''  fading.  Such 
fading  is  called  "selective"  because  it  occurs  only  for  frequencies  for  which  the  interference  condition  is 
satisfied.  At  other  nearby  receiver  locations  or  frequencies  the  time  delays  will  be  such  as  to  produce 
two  pulses  which  reenforce  each  other.  The  propagation  channel  can  be  simulated  by  a  multipath 
model  in  which  hT(t)  is  the  transmitted  wave  and  the  received  wave  is 

hR  it)  -  hT  it  -  t,),  (35) 

where  r ,  is  the  propagation  time  for  the  Ah  mode  or  hop  and  A,  is  the  amplitude  for  the  Ah  mode  or 
hop.  In  the  case  of  earth-ionosphere  waveguide  propagation  sometimes  the  waves  have  distinctly 
different  group  and  phase  velocities.  Then,  we  must  consider  both  the  phase  and  group  delay  times  r Xl) 
and  rXl).  The  phase  delay  of  the  baseband  modulation  is  governed  by  r p  \  and  the  group  delay  is 
governed  by 

THE  MSK  CORRELATION  RECEIVER 

One  can  readily  evaluate  the  performance  of  a  correlation  receiver  for  MSK  modulation  and 
coherent  detection  for  signal  reception  on  a  single-path  propagation  channel.  The  correlation  vector 
Rtirr)  obtained  when  one  uses  MSK  modulation  is  explicitly  calculated  by  substituting  the  expression 
for  hgiit)  from  Eq.  (32)  into  Eqs.  (27)  and  (28).  The  integral  /^(t,),  defined  as 

/*(t,)  -  Tjx  /frTy2  Aflj(r  +  t r)hg\it)  dt 
K- 1 

-  VI  U.  CZ'<’-*T>,U+'.-'T)  *  06) 

*— 0 
n'-O 

is  readily  evaluated.  The  limits  of  integration  in  this  equation  expend  far  beyond  the  region  of  variable 
t,  for  which  the  integrand  is  nonzero.  We  cair  change  the  integration  variable  to 

t'~  t-  nT.  (37) 

The  limits  of  integration  are  thereby  changed  to  tt  +  TJ1  -  nT  and  t,  -  Tg/2  -  nT. 

£-i 

/s<r,>  -  Til  I  /„/„  pit') pit'  +nT  +  r,-  nT)  dt.  (38) 
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Since  pit' )  —  0,  for  \t'\  >  T/2,  we  need  integrate  only  between  -  T/2  <  t  <  r/2.  We  obtain 

2  -+T/2 

lab,)  “  77*  I  /,/,  J.r/J  *(/')/>[/'  +  r,  +  in  -  n')T ]  dr.  (39) 


»-o 

»’-0 


(40) 


From  this  formulation  we  see  that  the  integral  depends  only  on  the  difference  between  n  and  n'and  not 
on  the  individual  values  of  n  and  n'  separately.  Let  us  define 

/T/2 

_T/2  pit')pit'  +  r,  +  kT)  dt. 

Since  k  -  n  -  n\  Eq.  (39)  becomes 

W  -  Tb~'  £  /,//,-,(r,) 

n—Q 
n'—0 


N  ...  N 
y-l  *+l-y 


-  T£' 

I 

I 

l, 

/*— 0 

k- n 

4-1 

N 

_  | 

7  1 

2 

~  T£' 

I 

4 

c 

«-0 

*—  1 

where 


/b/„-*Q(t,)  £/„_*. 


Un.k-l,  if0<  n-  k  <  1, 


)V 


-0,  if*— *<0orn  —  fc>— -1. 


We  then  interchange  the  order  of  summation  in  Eq.  (41)  to  get 


W  -  Tb 


-i 


C*(r,)  I  Inln-kU'-k. 

ff—0 


(41) 


(42) 


Because  of  the  orthogonality  of  the  set  { /„ }  with  any  sequentially  shifted  copy  of  itself,  we  can  write 


?- 


N 


I  Un-kUn-k-  f  8 


*,0. 


(43) 


«-o 


This  assumes  that  the  binary  sequence  {/„)  has  been  generated  by  a  method  which  has  "good*  auto¬ 
correlation  properties  for  the  pseudo-noise-coding  application  [7],  We  obtain 


NT~  * 

W  “  ~J~  C0(T,). 


To  calculate  Co(r,).  we  substitute  into  Eq.  (40) 

r  T/2 

C0irp)  “  J_r/,  pit') pit'  +  Tp)  dt- 


(44) 


(45) 


The  value  of  C0(rf)  is  symmetric  with  respect  to  rp  and  can  be  obtained  by  using  elementary  integra¬ 
tion  methods: 
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a  e  T/2-t 

C0It„)-  —  J  T  p(t')pW  +  rp)  dt.  if  tp  >  0. 

dt' 


2 

.  772— 


f  T 
2 


ir  t' 

cos 

t'  +  Tt 

T 

IT 

T 

_4 

T 


i  _  Jlii 

T 


irTp  ^  l  ■  Itp  I 

cos  — =-  +  —  sin  —— 
T  it  T 


Equations  (27)  and  (28)  become 


and 


In  matrix  form  we  have 


Xi(tp)  -  ( X,  cos  u>qt„  +  Y,  sin  <o0 t„)NTb  'C0(tp) 


Y,(t  p)  -  ( -X i  sin  <o0rp  +  Y,  cos  <o0t „)NTb  'Co(t,,). 


*>.)-  NTb'T' 


cos  o)qt p  sin  ai0Tp 

-sin&»oTp  cos  (Dor,, 


R,  C0( r  J. 


(46) 

(47a) 

(47b) 

(47c) 


A  receiver  can  be  adaptively  synchronized  to  the  incoming  signal  by  introducing  a  variable  delay 
tr  into  the  stream  of  baseband  keying  elements  to  generate  a  family  of  ideal  baseband  keying  elements 
hB i  and  hB1.  Mathematically  this  simply  means  the  introduction  of  the  variable  t  +  tb  into  the  expres¬ 
sion  for  hB  i  in  Eqs.  (27)  and  (28)  in  place  of  t.  The  selection  of  the  synchronization  time  tr  can 
greatly  affect  receiver  performance  when  the  actual  propagation  delay  time  t„  is  unknown  (or 
equivalently  the  clocks  that  generate  the  chip  streams  are  not  in  perfect  initial  synchronism).  Therefore 
we  rewrite  Eq.  (47c)  to  exhibit  the  dependence  of  system  performance  on  tr.  We  obtain  the  follow¬ 
ing  expression  for  R ,  (t,,t,): 


Ri(t„,tr) 


NTBy 


cos  (U0TP 
-sin  u>0Tp 


sin  m0Tp 
COS  ft»oTp 


Ri  Co(T p  ~  tr )• 


(48) 


The  magnitude  of  the  components  of  R,(,tp,tr )  are  proportional  to  the  in-phase  and  quadrature  com¬ 
ponents  of  the  received  radio-frequency  carrier  wave  because  of  the  multipliers  cos  w0tp  and  sin  <d0t„ 
in  the  matrix.  However,  the  magnitude  of  Rt(.Tp,  tr)  depends  on  the  propagation  delay  time  rp  and  on 
the  receiver  time  t*  by  the  factor  rp  -  tr.  If  the  receiver  is  badly  out  of  synchronization,  this  factor 
can  be  zero  or  very  small.  If  the  receiver  is  employing  phase-shift  keying,  any  change  in  the  propaga¬ 
tion  time  tp ,  even  after  the  receiver  is  perfectly  synchronized  on  the  chip  (  0  =  lrp  -  r^i  «  7),  will 
induce  a  corresponding  phase  change  in  the  phase  of  the  vector  Rt  -  (tp, tr)\  that  is,  the  doppler 
phase  variation  of  the  radio-frequency  carrier  induces  an  equal  doppler  shift  on  the  phase  angle  <t>R. 
These  types  of  doppler-signal  phase  variations  are  often  seen  when  an  mobile  transmitter  or  receiver  is 
used.  In  such  a  circumstance  we  cannot  assume  rp(t)  will  be  constant  at  each  instant  within  a  bit 
period  to  permit  the  factorization  shown  in  Eqs.  (21)  and  (22).  However,  the  other  manipulations  that 
we  performed  on  the  expression  for  R  are  all  valid  up  to  this  point.  To  obtain  the  resultant  correlation 
vector  Ri(tp,tr)  if  the  delay  time  is  slowly  changing  within  a  bit  period,  we  must  average  the  instan¬ 
taneous  values  of  RiWp(t),Tp]  over  an  entire  bit  period: 


S(T  r  >  -  NT-'T'  *T,n  costuo *,(') 
RArp,rR)~NTB  T  _sin<UoTf(,) 


S^lrt)  C»M'»  -  '»!• 


(49) 


From  this  expression  one  can  approximately  evaluate  the  affect  on  the  Rt  —  (rf(r).rK]  vector  for  vari 
ous  conditions  of  rpU)  and  synchronization  conditions  for  a  single  path  propagation. 
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We  have  reviewed  the  elementary  cases  of  a  stationary  and  moving  transmitter  having  a  single 
path  of  propagation,  because  this  is  the  simplest  case.  Linear  superposition  holds  for  the  waves  and 
correlation  procedures;  so  a  resultant  from  a  case  having  several  propagation  paths  can  be  obtained  by 
addition  of  terms  from  single-path  cases. 

TWO  PROPAGATION  MODES  AND  A  STATIONARY  TRANSMITTER 

When  a  stationary  transmitter  produces  a  signal  at  a  receiver  that  comes  by  two  propagation  paths, 
there  are  two  propagation  delay  times  r,,  and  t  p2  and  two  (real-valued)  amplitudes  of  the  two  modes 
A  |  and  A2.  Alternatively,  we  may  say  that  each  mode  is  described  by  a  complex  amplitude,  with  the 
magnitude  of  each  complex  amplitude  being  the  real  amplitude  and  the  phase  being  related  to  the  pro¬ 
pagation  time.  The  voltage  impressed  on  the  receiving  antenna  by  the  received  signal  represents  a 
waveform  that  is  the  sum  of  two  time-delayed  replicas  of  the  initial  waveform.  Let  the  subscript  i 
denote  the  mode  index.  Then  we  will  get  the  following  expression  for  the  total  correlation  vector 
RT(  tk): 

<- 1 


COS  U>0Tp  i 
-sin  <o0tp  , 


sin  ut0Tpi 

COS  O)0Tp  , 


Y,  co(tp,/  ~  T/f)- 


The  behavior  of  RI(tr)  as  a  function  of  the  receiver  chip  synchronization  time  tr  can  be  quite 
different  from  the  behavior  when  there  is  only  a  single  path.  A  splitting  of  the  correlation  peak  can 
occur.  Indeed,  instead  of  the  correlation  being  best  when  rp,-  tr,  it  could  be  bad  at  this  time. 
Numerical  examples  given  later  exhibit  this  circumstance. 

FIELD-STRENGTH  EXPRESSIONS 


For  the  case  of  an  antenna  trailing  behind  an  airborne  transmitting  facility,  the  correlation  vectors 
generated  by  both  the  horizontal  and  vertical  components  of  the  antenna  are  included  in  R?(tr).  For 
an  inclined  short  dipole  antenna,  the  endfire-generated  components  are  modified  by  cos  £  sin  y,  and  the 
broadside  generation  depends  on  sin  £  sin  y,  where  £  is  the  angle  between  the  horizontal  portion  of  the 
inclined  dipole  and  the  great  circle  path  to  the  receiver  and  y  is  the  inclination  angle  of  the  antenna 
measured  from  the  vertical.  The  expression  for  the  resultant  vertical  electric  field  Er(D)  at  distance  D 
from  the  transmitter  can  be  written  as  follows  (4,  p.  332J: 

Er(D)  -  cos  y  £"  +  sin  y  £*.  (51) 

In  this  expression 


—71  Ids 
h  VflAsin# 


£  S,' 5  A '  Gp(z,)Gf(z  )exp(ik0D  S„  -  in/ 4) 


—7)  Ids 
y/a  Asin0 


X  [a;  g;(z,)*;  (z.)  cos  £ 


rjA.*  Gp{z.)  ... 

+  — V'  sin<  exp  Uk0DS ,  -  i'ir/4). 

^ eh 
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The  quantities  in  these  expressions  are  explained  in  Ref.  4  and  are  summarized  below: 


A.' 

is  the 

A* 

is  the 

G'(z) 

is  the 

G*(z) 

is  the 

5, 

is  the 

z 

is  the 

2, 

is  the 

a 

is  the 

D  -  a9 

is  the 

A'(z) 

is  the 

A  *(z) 

is  the 

zfh 

is  the 

V 

is  the 

*  0 

is  the 

h 

is  the 

A 

is  the 

Ids 

is  the 

excitation  factor  of  the  qth  mode  for  a  vertically  polarized  electric  dipole, 
excitation  factor  of  the  gth  mode  for  a  horizontally  polarized  electric  dipole, 
height-gain  function  of  the  quasi-TM  vertical  electric-field  component, 
height-gain  function  of  the  quasi-TE  vertical  magnetic-field  component 
sine  i  '  the  complex  eigenangle  of  the  9th  waveguide  mode, 
height  of  the  receiver, 
height  of  the  transmitter, 
radius  of  the  earth, 

great-circle  distance  between  the  transmitter  and  the  receiver, 
normalized  impedance  of  the  quasi-TM-mode  fields, 
normalized  impedance  of  the  quasi-TE-mode  fields, 
coupling  impedance  in  the  mode  -  £#(0)///*(0)  at  the  earth’s  surface, 
impedance  of  free  space; 

free-space  propagation  constant  (-cu/c)  for  an  electromagnetic  wave  of  frequency 

reference  height  of  the  ionosphere, 
space  wavelength,  and 
current  moment. 


The  dispersion  that  occurs  on  each  individual  waveguide  mode  is  usually  minimal  when  well- 
behaved  types  of  D-  and  E-region  boundary  conditions  are  present  [61.  Each  individual  waveguide 
mode  can  be  considered  as  a  dispersionless  path  over  a  sufficiently  narrow  fractional  bandwidth.  The 
effect  of  propagation  is  to  produce  several  time-delayed  images  of  the  transmitted  waveform  at  the 
receiver.  The  amplitudes  of  these  waveform  images  vary  independently  with  time,  as  the  inclined 
antenna  is  transported  about  by  the  aircraft.  This  variation  with  position  of  the  aircraft  can  be  modeled 
readily  by  letting  the  amplitudes  M,)  and  delay  times  [tpJ  vary  in  an  appropriate  way. 


CORRELATION  EFFECTS  OF  ANTENNA  ROTATION 


It  is  easy  to  show  the  effects  of  rotation  by  an  inclined  elevated  antenna  given  the  previous 
development  of  the  mathematics.  The  horizontal  component  of  the  inclined  antenna  generates  both  TE 
and  TM  modes  according  to  Eq.  (S3),  because  the  horizontally  generated  fields  have  phases  and  ampli¬ 
tudes  different  from  those  generated  by  the  vertical  antenna  components.  The  waveguide  modal 
interference  pattern  for  an  elevated  horizontal  antenna  differs  from  that  of  a  vertical  antenna.  This 
effect  is  illustrated  in  Fig.  6  (taken  from  Ref.  3),  which  shows  the  fields  generated  by  variously  inclined 
antennas  at  20  kHz.  The  amplitudes  of  the  horizontally  generated  fields  vary  as  the  cosine  of  the  angle 
between  the  horizontal  antenna  component  and  the  direction  of  propagation. 

In  simulating  the  effect  of  this  variation  on  the  receiver  output  vector,  we  use  Eq.  (SO)  for  the 
received  correlation  vectors,  and  we  use  four  amplitudes  {A,),  four  phase  time  delays  {rf},  and  four 
group  time  delays  (rf).  Two  of  the  amplitudes  represent  two  vertically  generated  mode  fields  and  vary 
as  sin  y.  The  other  two  represent  horizontally  generated  TM  mode  fields  and  vary  as  the  product  cos 
{  cos  y.  The  resultant  expression  for  the  correlator  output  is  the  same  as  Eq.  (SO)  except  that  in  this 
case  the  correlation  output  is  the  summation  of  the  outputs  for  the  four  individual  modes. 

II 


*•  v. 
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Fig  6  —  Predicted  variation  of  field  strength  and  phase  with 
distance  for  20  kHz  and  vertical,  inclined,  and  horizontal  an¬ 
tenna  inclinations 


Denoting  the  X,(l>(rR)  from  the  vertically  electrically  generated  modes  as  X,<i>'"(tr)  and  the 
Xi(t>  (t  R  )  from  the  horizontally  electrically  generated  modes  as  T/')v(t„),  we  can  write 

N’  N1 

Xf(rR)  -  cos  y  £  */'>"  (r„)  +  cos  £  sin  y  £  X,M*  (r*).  (54) 

Parallel  definitions  can  be  made  for  Y^(tr),  K/')v,(tr)  and  YiU)i*(tr),  and  we  obtain 

Yf(rR)  -  cos  y  £  y/')v'(rg)  +  cos  £  sin  y  £  Y,u)hf  (r*).  (55) 

The  vector  Rj  (t*)  is  a  sum  of  two  vectors  Rj*  (rR)  and  Rj"  (t*): 

RI(,tr)  “  cos  y  R/'"  (r*)  +  cos  £  sin  y  R!*(tr).  (56) 

The  vectors  (t*)  and  Rl*  (rR)  can  be  parallel,  perpendicular,  or  opposed  to  each  other,  depending 
on  circumstances  and  the  value  of  tk.  The  resultant  R^  Wr)  vector  can  change  quite  dramatically  dur¬ 
ing  a  complete  orbit.  For  example,  if  Rjr  (rR)  and  R*  (rR)  are  of  the  same  magnitude  and_phase, 
there  would  be  one  large  oscillation  of  |/?;r(rR)|  per  orbit  (provided  ^  «  45  T  If  R"(tr)  and  RI*(tr) 
are  perpendicular  to  each  other,  there  would  be  two  oscillations  of  |/?r(T^)|  with  one  cycle  £.  There  is 
a  mathematical  isomorphism  between  the  phasor  E,(D)  and  the  vector  RJ{tr),  since  the  composition 
of  various  phasors  which  form  the  complex  phasor  sum  E,(D)  is  isomorphic  to  the_composition  of  vari¬ 
ous  correlation  vectors  which  form  the  two-dimensional  correlation  vector  sum  R^(tr).  When  both 
TE  and  TM  waveguide  mode  propagation  from  an  airborne  transmitting  facility  are  taken  into  amount, 
we  obtain 

R*(t r)  -  Rr(r  g)  cosy  +  sin y  lcos£  /?/*( TM, r  R )  +  sin  £ Rf*( TE,  t R ) ] .  (57) 
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The  total  vector  rJ(t  h  )  is  the  sum  of  independently  propagated  vertical  endfire-  and  broadside- 
generated  components.  These  components  are  themselves  sums  of  components  which  are  propagated 
by  different  waveguide  modes. 


DIAGRAMMATIC  REPRESENTATIONS 

It  is  instructive  to  view  diagrams  that  illustrate  the  behavior  of  the  various  Ri  vectors  as  r  K  and  ( 
vary.  The  following  example  illustrates  what  can  happen  in  a  particularly  severe  case  of  modal  interfer¬ 
ence.  Figure  7  depicts  the  variation  of  a  single-mode  correlation^ output  R/1)v,( tk)  from  a  stationary 
vertical  electric  transmitter.  Figure  8  shows  an  output  function  R/2,v*(t«)  for  a  signal  coming  slightly 
later  on  a  second-order  mode  and  180s  out  of  phase  with  the  first.  Figure  9  gives  the  resultant 
Rf'b  K),  which  is  the  sum  of  the  correlation  outputs  of  the  first  and  second  modes.  Figure  10  shows 
the  phase  of  the  vector  Rf'ir  K )  as  a  function  of  r  A. 


Fig  7  —  Graph  of  £/,)w<t*)  versus  r*  for  a  stationary 
vertical  electric-dipole  transmitter 


Fig.  8  —  Graph  of  /t,(J,w(Tr)  versus  r  K  for  a  stationary  vertical  electric  dipole 
transmitter.  The  phase  of  this  wave  is  180°  opposite  to  that  illustrated  in  Fig.  7. 
The  maximum  value  of  I  A/,,",(rg) I  occurs  at  slightly  greater  rx  values  than 
does  the  maximum  value  of  I  shown  in  Fig.  7. 
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Fit  10  —  Griph  of  the  antle  of  the  £,"(»•*)  vector 

with  the  x  axis  in  Fit,  9,  showint  an  abrupt  shift  of  phase  at 
r*-  0 
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Figure  11  shows  the  resultant  R^(t R)  —  cos  y  Rf'i rR)  +  sin  y  cos  {^/"(t*),  where  {  is  the 
angle  between  the  direction  of  propagation  and  the  endfire  direction  of  the  antenna.  The  angle  y  is  the 
inclination  angle  of  the  antenna;  we  take  y  “  45s  and  {  —  0  for  this  example.  Figure  12  shows  the 
same  resultant  except  for  the  condition j;  —  180°.  The  conditions  for  {  —  90°  and  {  —  270°  are  both 
shown  by  the  pattern  of  Fig.  9,  since  Rf{ t g)  is  composed  solely  of  the  vertically  generated  modes  in 
these  cases.  These  figures  give  a  picture  of  RI(tr)  at  the  four  cardinal  points  of  an  orbit. 


Y 


Fig.  II  —  Graph  of  S*(tr)  for  an  endfire  case  (cos  {  -  +1),  illus¬ 
trating  the  change  in  the  value  of  t R  for  which  I  |  -  0  and 

the  variation  in  the  size  of  the  peak  values  of  I  Kfh  g)  I  on  opposite 
sides  of  t„  —  0.  Here  a  relatively  small  value  of  I  R*h K)  I  is  as¬ 
sumed. 


Fig.  12  —  Graph  of  S,r(r  „)  fbr  the  opposite  endfire  case  (cos  {  -  -1)  to  the 
case  in  Fig.  II.  Here  the  largest  value  of  I  R/r(r k )  1  occurs  for  r„  >  0 
whereas,  in  Fig.  II  it  was  for  t*  <  0.  Synchronization  times  (  rN)  to  be  con¬ 
sidered  in  subsequent  figures  are  shown. 
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In  Fig.  12  four  sample  synchronization  times,  designated  t#,.  tR2,  rA],  and  t*4,  represent  possi¬ 
ble  correlation  times  that  might  occur  within  a  receiver.  We  look  at  the  variation  of  the  correlation 
vector  magnitude  IaF/Vt*,)!  and  phase JiJ(tRi)  during  a  complete  orbit  for  these  times.  Figures  13  and 
14  show  the  phase  and  magnitude  of  Ri(rKt).  The  phase  is  constant  in  this  case,  because  is 

greater  than  and  parallel  to  The  magnitude  varies  in  a  mild  way.  In  Figs.  IS  and  16  showing 

the  phase  and  amplitude  for  the  time  rX2,  the  phase  makes  two  quick  successive  180°  changes  in  the 
neighborhood  of  (  -  180°,  and  the  amplitude  goes  through  two  successive  nulls.  Figures  17  and  18 
give  R,r( 0).  This  vector  varies  cosinusoidally,  because  there  is  no  contribution  from  the  vertical  elec¬ 
tric  antenna  component.  Figures  19  and  20  show  the  behavior  of  RI(trj),  which  is  similar  to  the 
behavior  of  /?/r(r*2)  except  that  the  nulls  occur  near  £  —  0°  rather  than  £  -*  180°.  Figures  21  and  22 
show  the  behavior  of  Rf(rR4),  which  is  similar  to  except  for  a  shift  in  orbit  angle  £  of  180°. 


=  tR1) 


Fig  13  —  Graph  of  4/r(rxl)  as  the  angle  of  turn  {  of  an  inclined  antenna 
changes  from  —180°  to  +  180'  for  the  synchronization  time  tr[.  This  phase 
is  constant  because  the  contribution  of  the  vertical  antenna  component  always 
exceeds  that  of  the  horizontal  component  for  this  synchronization  time. 


Fig.  14  -  Graph  of  l#,r(»>|)|  as  ( 
varies  from  —180*  to  +180'.  The  rela¬ 
tively  gentle  variation  of  lff,r(r4|)l  oc¬ 
curs  because  I  I  >  I  | )  I. 
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Fig.  15  —  Graph  of  ^/r(r#2)  versus  orbit  angle  {.  Fo  Fig.  16  —  Graph  of  I  Rj(r  k 2)  I  versus  orbit  angle  {,  show- 

this  lock-on  time  there  are  two  rapid  and  large  phase  ing  how  this  function  goes  to  zero  when  42r(r*2)  is  making 

changes  near  l  —  180°  the  abrupt  phase  transition 


Fig.  17  —  Graph  of  -  0)  versus  angle  {. 

There  is  essentially  no  contribution  from  the  vertical 
component  for  this  value  of  tk,  so  the  phase  is  con¬ 
trolled  by  the  horizontal  component  of  the  transmit¬ 
ting  antenna. 


Fig.  I*  —  Graph  of  I  R/(t*  -  0)|  versus  orbit  angle  £ 
The  cosinusoidal  variation  is  caused  by  the  cos  £ 
modulation  of  the  horizontal  component.  The  vertical 
component  contributes  essentially  zero. 


Fig.  19  —  Graph  of  */r(r#j)  versus  orbit  angle  (.  For  this  syn¬ 
chronization  time  the  abrupt  changes  in  phase  occur  around  the 
(  -  0  point  rather  than  around  the  {  «  180*  for  the  r„2  case  il¬ 
lustrated  in  Fig.  1 5. 


ig.  20  -  Graph  of  lR/r<r*j  ^versus  orbit  an- 
e,(.  The  smallest  values  of  lR/r(r„2ll  occur  for 
rbst  angles  near  zero  degrees  rather  than  near 
JO*  as  was  the  case  for  l#,r<r*2)l  »*  shown  in 
ig  10 
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Fig.  21  -  Graph  of4,r(TM4)  versus  orbit  angle  {. 
For  this  receiver  time  the  vertically^  generated 
correlation  vector  dominates;  that  is  > 

I  No  change  of  phase  can  occur,  since 

Kf  .-*4)  and  are  parallel. 


Fig  22  —  Graph  of  lS,r<x  M)1  versus  orbit  angle 
{.  For  this  receiver  time  the  maximum  correla¬ 
tion  vector  occurs  for  {  -  180",  whereas  for  t  „ , 
the  maximum  correlation-vector  magnitude  oc¬ 
curred  at  {  -  0°. 

Finally,  Fig.  23  shows  an  estimate  of  rRofl,  the  optimum  value  of  receiver  correlation  time  for 
different  values  of  orbit  angle  £.  This  graph  shows  that  r*,  and  rff4  are  preferable  when  |/Mrjti)l  and 
|£,(tr4)|  are  largest.  The  phase  is  constant,  and  the  amplitude  varies  gently.  Relatively  mild  ampli¬ 
tude  behavior  occurs  for  this  correlation  time  because^ I R/v'(ts1)|  >  |R,**(tjh)|.  The  phase  js  constant 
because  /F,y,(r*,)  is  parallel  to  Although  r!"(tri)  can  change  the  magnitude  of  Rf(rRl),  it 

does  not  change  the  phase. 


Fig.  23  —  Graph  showing  ihe  values  of  tr  that  should  be 
chosen  for  various  orbit  angles  (  for  the  case  illustrated  in 
Figs.  7  through  22 
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In  the  preceding  example,  there  was  minimum  correlation  near  the  middle  of  the  correlation 
envelope,  and  the  correlation  maximum  was  greatly  displaced  from  the  maximum  point  of  the  correla¬ 
tion  of  either  individual  multipath  signal.  Thus,  the  rather  small  difference  in  arrival  times  r,,  -  r,2 
changes  the  optimum  chip  synchronization  time  tk  by  a  major  amount.  This  effect  is  fairly  easy  to 
explain.  The  condition  described  by  our  equations  is  a  severe  case  of  multipath  interference.  Two  sig¬ 
nals  come  to  the  receiver  antenna  offset  in  time  by  r,,  -  r,2.  The  first  and  second  signals  travel 
independently  of  each  other,  so  when  the  first  signal  pulse  arrives,  the  signal  begins  to  build  up.  Then 
the  second  signal  arrives  out  of  phase  and  cancels  out  the  first  signal  for  a  long  interval  of  time  while 
both  are  strong.  Then  the  first  signal  goes  to  zero  because  its  waveform  has  passed  by,  while  the 
second  signal  remains  to  give  an  excitation  of  the  antenna  again.  Thus  the  strongest  portions  of  each 
signaling  element  in  this  case  of  interference  occur  at  the  beginning  and  the  end  of  each  signaling  ele¬ 
ment.  In  the  middle  of  the  element,  the  two  signals  arriving  from  different  paths  interfere,  and  there  is 
a  reduction  of  the  field  present. 

When  a  receiver  tries  to  synchronize  with  the  signal  at  its  strongest  point,  it  must  search  for  syn¬ 
chronization  either  at  the  beginning  of  the  chip  or  at  the  end  but  not  in  the  middle.  Thus,  if  a  system 
were  presynchronized  with  an  optimum  r  R  for  either  of  the  individual  mode  delays  or  for  an  average  of 
the  two  mode  delays,  it  would  actually  be  poorly  synchronized  for  optimum  operation!  For  better 
operation  three  correlators  might  be  used— one  adjusted  at  the  central  expected  time  of  arrival  of  the 
wave,  and  one  at  each  end  of  the  chip.  If  the  two  waveforms  reinforce  each  other  to  give  the  max¬ 
imum  in  the  middle,  the  central  correlation  peak  will  be  large  and  the  two  side  peaks  will  be  small.  If 
however,  an  interference  condition  is  present,  the  middle  correlator  will  be  weak  and  the  two  side 
correlators  will  be  strong. 

NUMERICAL  CALCULATIONS 

Procedures 

We  will  calculate  some  values  of  Rj  for  cases  that  are  inspired  by  realistic  mode-theory  calcula¬ 
tions  of  the  phase  and  amplitude  of  VLF  waves.  For  the  calculated  fields  shown  in  Fig.  6,  the 
waveguide  mode  constants  are  given  in  Table  1 . 


Table  1  —  Constants  Used  in  the  Calculations 


A,  -  0.4 

G,(z)-  1 

A2-  1.2 

G2(z)  -  1 

S i  -  1.002  +  0.000689i 

A,(z,)  -  -0.02i 

Si  —  0.003  +  0.001 38i 

Mz,)-  0.05i 

7-  0  km 

r0  -  4  km 

G,(z,)-  0.965  +  0.259i 

z,  -  10  km 

Gi(zt)  -  0.95 

Radiated  power  —  1  kw 

These  constants  were  used  in  the  following  equations  to  generate  the  field  -strength  values  in  Fig.  6: 


-  iWe0)l/2/ 

(57) 

and 

El"(r,0)  -  /(^0/e)l/2/  <*vcos$  ]£  A„(z,)S®5F,. 

(58) 

where 

f *  (da  j,) 1/2  ( sin  (D/a ) )  G" (z)  G" (z,)A"  exp  (,ir/4  +  lk°  DS") 

(59) 
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These  equations  are  simplified  versions  of  Eqs.  (5U  through  (53),  since  they  contain  only  TM- 
mode  constants  of  an  isotropic  ionosphere.  We  calculate  Rlil)ve  for  each  mode  according  to 


x,u)ye  (t*,£)  -  N\Ar\  cos(to0Tf-'"(£)l  C  (t r‘(0  -  rR) 

(60) 

and 

(tr.O  -  NUr\  sin lw0rf' *'(£))  C  (rf  «(£)  -  t,). 

(61) 

where 

N  is  a  normalizing  constant, 

is  amplitude  of  the  Ah  TM  mode  generated  by  the  vertical  component  of  the  transmit¬ 
ting  antenna, 


«»o  is  the  angular  frequency  of  the  wave, 

r /’•'"({)  is  the  delay  time  of  the  phase  of  the  wave  when  the  aircraft  is  at  angle  £  of  its  orbit  for 

the  Ah  waveguide  mode,  and 

Tf  vr(()  is  the  group  delay  time  of  the  wave  when  the  aircraft  is  at  angle  £  of  its  orbit  for  the  Ah 
waveguide  mode. 

The  amplitude  is  given  by 

A?‘-  S,xiG,(z)  Gl(zs)\,  exp (iir/r  +  i  k0  DSf).  (62) 

The  phase  delay  time  rf  vr(£)  is  obtained  by  evaluating  the  phase  of  A,: 

a>orre(0-  argCVO.  (63) 

The  group  delay  time  r  f- v,(£ )  is  obtained  by  evaluating  the  derivative  of  the  phase  of  A,\ 

Tf  v'(£)  -  -^largCV')].  (64) 

The  path  length  D  from  the  transmitting  antenna  to  the  receiver  location  is  a  function  of  the  orbit  angle 
£  and  the  orbit  radius  r  according  to  the  approximation 

D~D0+r  sin  £,  (65) 

where  D0  is  the  distance  from  the  receiver  to  the  center  of  the  orbit.  To  calculate  /?"(£, r*),  we  per¬ 
formed  summations  as  indicated  in  Eqs.  (54)  through  (56). 

We  calculate  the  horizontally  electrically  generated  fields  in  an  analogous  manner: 

X,i,)ht  (r*.£)  -  N\Al"\  coslw0Tf  v(£)l  C0  (rf  *(£)  -  t  R)  (66) 

and 

y,U)he  (r*.£)  -  NUf  I  sin {<u0T /’•*'(£)]  C0  (rf*(£)  -  r*).  (67) 

The  amplitude  of  the  i  th  mode  is  given  as  follows: 

A*  -  G,(z)G,(z,) A,  exp(nr/4  +  /*„  D  S,)  S, 5  A,(z,)  -  A?'  St  A ,(z0)-  (68) 

The  phase  and  group  delay  times  are  given  in  the  same  manner  as  for  the  vertically  generated  modes: 

«0rf*<£)-  argM/*)  (69) 
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and 

rf*<£)-  ^-larg04,*0].  (64) 

00) 

Tables  2  and  3  give  the  numerical  values  used  in  making  the  correlation  calculations. 


Table  2  —  Set  of  Quantities  Used  in  Calculation 
of  Correlation  Vectors 


Mode 

Quantity 

Value 

1 

MH 

Tf*' 

Tf« 

1  A?\ 
rf* 
rf* 

1.26500262  x  10” 1 
1.33658211  x  ir2 
1.34400  x  10r2 
2.52495471  x  1<T3 
1.33533157  x  ir2 
1.34400  x  ir2 

2 

Mri 

rf« 

Tf" 

Mf  1 

Tf* 

Tf* 

1.11773228  x  10"1 
1.32437617  x  10'2 
1.356666  x  10~J 
5.62805237  x  10~3 
1.32562578  x  ir2 
1.356666  x  ir2 

Table  3  —  Second  Set  of  Quantities  Used  in 
Calculation  of  Correlation  Vectors 


Mode 

Quantity 

Value 

1 

Ur\ 

Tf" 

Tf" 

MM 

Tf.* 

Tf* 

1.26500262  x  10r> 
1.33658211  x  ir2 
1.34400  x  1(T2 
2.52495471  x  ir3 
1.3315820  x  ir2 
1.34400  x  ir2 

2 

Mr  i 

Tf" 

Tf" 

MM 

Tf  * 
Tf* 

1.26500262  x  ir1 
1.33908211  x  10-2 
1.356666  x  10~2 
5.62805237  x  ir3 
1.3365820  x  l<r2 
1.356666  x  ir2 

Results 

Figures  24  through  29  show  rt,r(r*)  and  0,r(7>)  contours  for  y  values  of  90°,  85°,  80°,  60°,  45°, 
and  0°  with  the  propagation  parameters  of  Table  2.  Figure  30  shows  #?(,)"(t*),  the  correlation  vector 
of  the  first-order  mode.  Figures  31  through  36  show  the  correlation  vector  for  propagation  parameters 
of  Table  3  for  y  values  of  90®,  85®,  80®,  60®,  45®,  and  0®.  The  "average*  chip  duration  for  each  of  these 
calculations  was  rf  -  1  ms.  The  total  length  of  each  chip  is  f  -  2tc  -  2  ms.  The  radius  of  the  orbit  is 
3  km;  the  distance  from  the  center  of  the  orbit  to  the  receiver  is  4000  km. 
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>a 


Fig.  24a  —  Amplitude  behavior  of  R*(t k )  for  y  —  90°, 
calculated  with  the  parameters  of  Table  2 


Fig.  24b  —  Phase  behavior  of  »>)  Tor  y  —  90°, 
calculated  with  the  parameters  of  Table  2 
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Fig.  25a  Amplitude  behavior  of  Rf(TK)  for  y  —  85° 
calculated  with  the  parameters  of  Table  2 


Fi*.  25b  -  Phase  behavior  of  ^r(r*)  for  y  -  85° 
calculated  with  the  parameters  of  Table  2 
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Tr(  rnULtytctmcb) 


Fig.  25d  —  Phase  behavior  corresponding  to  Fig.  25b 
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Fi*.  28*  —  Amplitude  behavior  of  for  y  —  45°, 

calculated  with  the  parameters  of  Table  2 
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R'  (Tr  ) 


Fi*.  28c  —  Amplitude  contours  corresponding  to  Fig.  28a 
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Fig.  29a  —  Amplitude  behavior  of  R,r(Tj)  for  y  m  0°, 
calculated  with  the  parameters  of  Table  2 


r 

Fig.  29b  -  Phase  behavior  of  for  y  -  0°. 

calculated  with  the  parameters  of  Table  2 
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Fig.  3 Id  -  Phase  behavior  corresponding  lo  Fig.  31b 
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—  Phase  behavior  of  M,t(tk)  for  y  *■ 
calculated  with  the  parameters  of  Table  3 
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Phase  behavior  corresponding  to  Fig.  32b 
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Fig  36«  —  Amplitude  behavior  of  £,r(i>)  for  y  -  0°, 
calculated  with  the  parameters  of  Table  2 
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Fi*.  36c  —  Amplitude  contours  corresponding  to  Fig.  36a 


10  II  12  13  14  t»  is 

Fig.  36d  —  Phase  behavior  corresponding  to  Fig.  36b 


KELLY,  HAUSER,  BECK.  AND  RHOADS 


Dlscmstaa 

Figure  24  shows  the  correlator  performance  for  the  case  of  a  completely  horizontal  antenna.  For 
this  case  I/?/(t*)|  goes  to  zero  at  £  —  90°  and  270°,  and  there  are  sudden  180°  phase  flips  at  these 
angles.  The  best  correlation  time  appears  to  be  at  tk  -  13.S  ms,  which  is  the  middle  of  the  correlation 
peak. 


Figure  25  shows  the  correlator  performance  for  the  case  of  an  antenna  with  inclination  y  **  85° 
from  the  vertical.  For  this  case  it  is  difficult  to  say  if  there  is  a  single  best  receiver  synchonization  time 
tr.  During  most  of  the  orbit  it  appearsjhat  the  time  tr  —  13. 5  ms  would  be  the  best.  However,  this 
value  of  t/j  appears  to  give  a  value  of  |/?,r(rR)|  that  is  severely  reduced  at  £  —  90°  and  £  —  270°.  For 
this  case  a  system  of  selection  of  three  values  of  tr  separated  by  1-ms  intervals  (12.5,  13.5,  and  14.5 
ms)  would  give  a  more  stable  average  performance  in  the  evaluation  of  Rf(TK).  The  discontinuity  in 
Fig.  25d  at  r*  —  13.8  ms,  signified  by  the  heavy  vertical  line,  is  due  to  a_360°  ambiguity  which  has  no 
real  significance  except  at  the  nulls  (£  —  120°  and  £  “  240°),  where  |/?/r(rR)|  — *  0.  Comparison  of 
Figs.  25c  and  25d  illustrate  this  feature  clearly. 

Figure  26  shows  the  correlation  performance  for  an  antenna  with  y  —  80s.  In  the  vicinity  of  t* 
—  14.0  ms  and  £  —  180°  there  is  an  abrupt  drop  of  signal  and  jump  in  phase.  Otherwise  the  |A/(ta)| 
function  seems  fairly  well  behaved. 

Figure  27  shows  the  behavior  of  Rfi t  k  )  for  an  antenna  with  y  ~  60s.  The  amplitude  and  phase 
variation  of  the  correlation  are  much  milder  than  for  the  cases  with  smaller  values  of  y,  since  the  verti¬ 
cal  electric-field  components  exceed  the  horizontal  components  by  it  least  one  order  of  magnitude  at 
this  droop  angle.  This  trend  is  continued  as  y  is  decreased  to  45s  and  then  to  0s  in  Figs.  28  and  29. 

Figure  30  gives  the  value  of  R/I)v'(r  „).  It  illustrates  the  case  of  pure  first-order-mode  propaga¬ 
tion  from  a  vertical  electric-dipole  source. 

The  values  of  the  propagation  quantities  in  Table  3  are  similar  to  those  of  Table  2  except  for  an 
adjustment  to  produce  maximum  mode  cancellation  and  interference  effects.  We  have  made  the  first- 
order  and  second-order  vertically  generated  modes  have  opposite  phase  and  equal  magnitude.  We  have 
modified  the  phase  of  the  horizontally  generated  modes  so  they  are  in  the  same  phase  plane  with  the 
vertically  generated  modes  to  produce  maximum  field-strength  variation.  Hence  both  severe  effects 
appear  in  Figs.  31  through  36. 

Figure  31  gives  the  behavior  of  for  the  propagation  parameters  of  Table  3.  The  results 

are  almost  identical  to  those  of  Fig.  24,  since  y  -  90s  and  the  only  changes  made  in  the  parameters  of 
the  horizontally  generated  modes  are  to  the  phase  delay  times  r  f  v  and  r  f v. 

Figure  32  shows  the  behavior  of  R*( rK)  for  the  case  of  y  »  85s.  The  influence  of  the  splitting 
of  the  vertically  generated  function  is  more  evident  on  the  phase  plot  than  on  the  amplitude  plot. 

Figure  33  shows  the  behavior  of  R*(tr)  for  the  case  of  y  “  80s.  Here,  the  influence  of  the 
correlation-peak  splitting  of  the  vertically  generated  function  is  quite  evident  as  the  notch  in  the  corre¬ 
lation  peak  makes  its  way  from  side  to  side  through  the  values  of  r  K  and  £ .  This  is  illustrated  further 
in  Figs.  34  and  35,  which  give  the  graphs  for  y  —  60°  and  45s. 

Finally,  Fig.  36  shows  the  split  correlation  peak  for  a  purely  vertical  antenna  being  carried  around 
an  orbit.  Here  we  find  a  nearly  constant  variation  of  |/?(r(r,)|  versus  £,  and  the  phase  of  the  correlation 
function  41(tr)  varies  only  moderately  with  £.  As  a  function  of  rR,  these  variables  show  the  expected 
splitting. 
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The  phase  changes  gradually  with  tr  in  Figs.  27  through  29,  whereas  in  Figs  34  through 

36  it  changes  abruptly  at  certain  tr  values  but  then  achieves  a  constant  value  again  on  either  side  of 
the  change  point.  This  is  readily  explained  using  the  diagrammatic  techniques  previously  discussed, 
since  in  the  first  case  the  correlation  function  for  the  individual  modes  /?/*’( tr)  are  not  coplanar,  and 
in  the  second  case  they  are  coplanar,  so  that  only  abrupt  phase  changes  are  permitted. 

CONCLUSIONS 

This  report  shows  possible  types  of  variation  of  the  correlation  vector  Rfir R)  when  an  airborne 
VLF  transmitting  system  is  employed.  Use  of  a  single  correlator-synchronization  time  can  result  in  loss 
of  communications  because  of  possible  signal  degradation  and  phase  change  caused  by  orbital  and 
modal  interference  effects.  Multiple  and/or  adaptive  correlation  and  decision  techniques  seem  indicated 
for  this  channel. 
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